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We examine the MiniBooNE neutrino, MiniBooNE antineutrino and LSND antineutrino data 

sets in a two-neutrino ¥~v^ e oscillation approximation subject to non-standard matter effects. 
We assume those effects can be parametrized by an //-independent effective potential, V s = ±A S , 
experienced only by an intermediate, non-weakly-interacting (sterile) neutrino state which we assume 
participates in the oscillation, where +/— corresponds to neutrino/antineutrino propagation. We 
discuss the mathematical framework in which such oscillations arise in detail, and derive the relevant 
oscillation probability as a function of the vacuum oscillation parameters Am 2 and sin 2 2# Me , and 
the matter effect parameter A s . We are able to successfully fit all three data sets, including the 
MiniBooNE low energy excess, with the following best-fit model parameters: Am 2 = 0.47 eV 2 , 
sin 2 20 Me = 0.010, and A s = 2.0 x 10~ 10 eV. The ^-probability for the best fit corresponds to 
21.6%, to be compared to 6.8% for a fit where A„ has been set to zero, corresponding to a (3+1) 
sterile neutrino oscillation model. We find that the compatibility between the three data sets 
corresponds to 17.4%, to be compared to 2.3% for A s — 0. Finally, given the fit results, we 
examine consequences for reactor, solar, and atmospheric oscillations. For this paper, the presented 
model is empirically driven, but the results obtained can be directly used to investigate various 
phenomenological interpretations such as non-standard matter effects. 



I. INTRODUCTION 



Neutrino flavor oscillation is a phenomenon that arises due to non-zero, distinct neutrino masses and lcptonic mixing. 
The current Standard Model picture incorporating neutrino oscillation relates three neutrino mass eigcnstatcs with 
masses mi, to 2 and m.3 to the three neutrino weak flavor eigcnstatcs, v e , v^, and v T . This is done through the lcptonic 
mixing matrix, U , defined as a product of three rotations, 
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Using the above definition, neutrino flavor change can be described as a function of the mixing elements and neutrino 
masses in terms of the three-neutrino oscillation probability 




P[y a ^up) = 6 a p ~ AY,^{U* m U aj U Pl U* P] }sm 2 (1.267 AmlL/E) 

i>j 

+ Y^^{KiU aj U0 i U^}sm(2.53AAm 2 ij L/E) , (2) 



where a, /? = e, /i, t; i, j = 1, 2, 3; Am,, = mf — mA in eV , E is the neutrino energy in MeV; and L is the neutrino 



propagation distance (in the lab frame) from production to detection in meters. In a simpler, two-neutrino oscillation 
approximation, the oscillation probability is reduced to 

P{ { v\-^v\) = S al3 - AlU^U^U^U^lsin 2 (1.267 A^L/E) . (3) 
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The above approximation holds when Ato 2 ^ and Am 2 ^ differ significantly, and when only two of the three weak flavors 
(a and 0) participate in the oscillation most-dominantly. 

Neutrino oscillations described by the oscillation probabilities in Eqs. [2] and [3] have now been established through 
multiple experiments [l[ ■ Those experiments study two- neutrino oscillations in the form of appearance of neutrinos 
of flavor /3 in a neutrino beam of flavor a, described by Eq. [3] when f3 ^ a (S a p = 0), or in the form of disappearance 
of a neutrino beam of flavor a, described by Eq. [3] when a = /3 (8 a p = 1). 

Almost all available experimental results are consistent with the following vacuum oscillation parameters [l[: 

Amij =m\-m\ = 7.65 x 1CT 5 cV 2 , 

sin 2 0i2 = 0.304 , 
Am 2 : = m 2 - to? = 2.40 x 10" 3 cV 2 , 
sin 2 #23 = 0.5 , 

< sin 2 26»i3 < 0.15 . (4) 

There are, however, existing results which cannot be accommodated within this picture, and suggest a possible need 
for extension beyond the three-neutrino framework. Those results come from both appearance and disappearance 
measurements performed at relatively short baselines. 

More specifically, two independent experiments, LSND and MiniBooNE [a-Q, have observed three independent 
appearance-like excesses of electron neutrinos and/or antineutrinos in muon neutrino and/or antineutrino beams, 
at least two of which are consistent with oscillations at the level of 2.8-3.8 a. Under a two-neutrino oscillation 
approximation, each of those two measurements reveal excesses which correspond to a large Ato 2 , such that Ato 2 3> 
Am| 2 S> Am^i, requiring at least one extra neutrino mass eigenstatc to be added to the standard three-neutrino mass 
spectrum. This extra mass eigenstate is assumed to be mostly "sterile" , i.e. non-weakly-interacting, by assumptions 
of unitarity and experimental constraints from Z — >• vv decay measurements Q. Such models are referred to as 3 
active + 1 sterile neutrino models, or, (3+1), and are usually explored in a two-neutrino oscillation approximation. 
In the case of MiniBooNE and/or LSND, this is done by setting Ato| 2 ~ Am|i — and Am|j ~ Am 2 3 ~ Am 2 2 , so 
that the appearance oscillation probability can be obtained from Eq. [3j for a = /x, /3 = e, and Am^ = Am 2 i, 

p(v\-^u \) = sin 2 20 Me sin 2 (1.267Am^L/£:) , (5) 

where we have used the orthogonality relation U^iU^ = —U^tJJ*^ and definition sin 2 28^ e = 4|£/ e 4| 2 |£/ Al 4| 2 . 

At the same time, a number of short-baseline reactor antineutrino experiments have looked for electron antineutrino 
disappearance at the same range of Am^ , described by the v e survival probability 

P{ [ v ] e -* v \) = I -sm 2 29 ee sm 2 (l.27Aml 1 L/E) 

= 1 - 4|t/ e4 | 2 (l - \U ei \ 2 ) sin 2 (1.27Am 2 iL/i?) , (6) 

where we have replaced |t4i| 2 with (1 — |Z7 e 4| 2 ), by way of unitarity. While those same short-baseline reactor disap- 
pearance searches have provided strong limits in the past, a recent re-analysis of predicted reactor antineutrino fluxes 
revealed an underestimation of flux predictions previously assumed by those experiments (Io| . As a result, the limits 
on sin 2 29 ee have now been evaded, and an overall normalization reduction observed in their spectra shows consistency 
with sin 2 26 ee ~ 0.1 and Am^ > 1 eV 2 [ll|. 

It is tempting to attribute the MiniBooNE, LSND, and reactor short-baseline signals to the existence of a single 
sterile neutrino; however, attempted (3+1) fits have demonstrated that at least MiniBooNE neutrino results and LSND 
antineutrino results are incompatible under this scenario [l2| , suggesting that any successful theoretical interpretation 
of those results must be more complex than simply (3+1). The source of this incompatibility is the fact that the 
MiniBooNE antineutrino mode and LSND antineutrino mode appearance searches prefer moderately small sin 2 29^ 
and Am 2 ~ 1 cV 2 , while the MiniBooNE neutrino mode appearance search disfavors such sin 2 26^ values for the 
same Am 2 ~ 1 eV 2 . 

A minimal extension to the (3+1) model, which would allow for CP-violation (i.e. P{v^ — > v e ) ^ P(p^ — > D e )), 
can successfully fit all of the above three signatures and seems plausible as an explanation. However, it still requires 
relatively large mixing amplitudes in order to reasonably accommodate them, and so it is consequently disfavored by 
disappearance experimental constraints p^l - [l5| . 

In view of the shortcoming of the above, or "vacuum," sterile neutrino oscillation models in reconciling MiniBooNE 
and LSND results, phenomenological efforts have now turned toward consideration of CPT- violating models [l6l - [2ll |. 
or effectively CPT- violating models which involve non-standard matter effects [2214251 ] . Motivated by the latter class 
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of models, in this paper, we consider a four- neutrino oscillation scenario, where the fourth neutrino flavor state, v s , is 
subject to matter effects due to some interaction potential of the form 

V s = ±A S , (7) 

where A s is a constant and +/— corresponds to neutrino/antincutrino propagation through matter. This gener- 
alization allows us to follow an agnostic approach as to the underlying source of this effect. We assume that the 
excesses observed by MiniBooNE and LSND are manifestations of — > v e oscillation via the fourth mass eigenstate, 
which is assumed to be on the order of 0.01-100 cV 2 , and a mix of v e and eigenstates at the < 5% level each 
(|C7 e4 | 2 , \U^\ 2 < 0.05), and v s at the > 90% level (|J7 s4 | 2 > 0.90). 

We discuss the oscillation framework in detail in the following section. In Sec. IIIII we describe the analysis and fit 
machinery used to apply this framework to the MiniBooNE and LSND data sets. In Sec. HVl we present quantitative 
results, first in a (3+1) scenario without this matter effect (^4 S = 0), as a reference, and then with the matter effect 
turned on. In Sec. IIV1 we provide a qualitative description of the results, and, in Sec. I VII we discuss implications for 
atmospheric, solar, and reactor experiments. Conclusions arc presented in Sec. IVII1 



II. OSCILLATION FRAMEWORK 



In this section, we derive the Vn — > v e and — > v e appearance oscillation probabilities to which we attribute the 
observed MiniBooNE and LSND excesses. For simplicity, we assume that the matter potential V s experienced by v s 
is much larger in amplitude than the standard model matter effect potentials, Vcc and Vnc, so that the effective 
matter potential in neutrino flavor space can be approximated as 
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The effective Hamiltonian for neutrino propagation in matter, expressed in flavor space, is given by 



H m = — U 
2E 






m\ 



+ V 



(9) 



For simplicity, we have assumed that m 2 , m\ and m| in vacuum are degenerate and negligible relative to m 2 , with 
AmW = Am 2 being the only mass squared difference in vacuum. We also assume that m 2 /2E and m\/2E are 
negligible relative to V s . Then, using the standard form of the U mixing matrix in vacuum, 
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the effective Hamiltonian in matter becomes 
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Diagonalizing the effective Hamiltonian in Eg. 1111 gives the following eigenvalues: 

Ai = , 
A 2 = , 



A 3 = ^ (2EV S + Am 2 £ Q \U a4 \ 2 - ^{-2EV S - Am 2 J2 a \U a4 \ 2 ) 2 - 8£1/ s Ato 2 (£ q |[/ q4 | 2 - \U si \ 2 )) 
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4E 
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\U S 4\ 2 )) 



(10) 
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which, by unitarity, arc reduced to 



Ai = , 
A 2 =0 , 



X 3 = je [ 2EV * + Am2 - V( 2EV s + Am 2 ) 2 - SEV s Am 2 {\ - |£/ s4 | 2 ) 
A 4 = (2EV S + Am 2 + y / {2EV s + Am 2 ) 2 - 8EV s Am 2 (l - |t/ s4 | 2 )) ■ (13) 

The differences A 4 — Ai,2, A 4 — A3, and A3 — Ai^ suggest three distinct effective Am 2 M values in matter. Note, however, 
that, if the active flavor content of the fourth mass eigenstate is small, then we can approximate 

l-|[/ s4 | 2 ~0, (14) 

in which case all eigenvalues except A 4 become zero. The non-zero A 4 , 

X A = ^(2EV S + Am 2 ) , (15) 

implies one effective Am 2 Ml which corresponds to 

Am 2 M = Am 2 + 2EV S . (16) 

We will be using the above Am 2 M approximation in our fits, which is a justified assum ptio n according to the level of 
unitarity in the three-neutrino mixing matrix, which is experimentally constrained |26rl28|. Note that, when V s —> 0, 
Am 2 M reduces to the vacuum Am 2 value, as expected. 

The matrix consisting of the unit-normalized eigenvectors (columns) of the effective Hamiltonian in Eq. 111! defines 
the new mixing matrix in matter, and respective mixing elements, 

U M = P . (17) 

For the purposes of this paper, we are interested in the general expression from which the — > v e and — > u e 
appearance oscillation probabilities are derived, 

P(v» -+ ve) = I £ U* el e-^ L ' 2E U, t \ 2 . (18) 



According to Eq. [TBI the elements of interest to this appearance channel are U^f and U^j . Since only one Am 2 

jM 
' /'4 



dominates, that associated with A 4 , the parameters of interest are just (f/Jf)* and UM. From the eigenvectors of H m , 



one identifies 



M 2 * / l+4Am 2 |a s4 ! 2 (l-|y s 4| 2 )/(2By s -ATn 2 (l-2|C/ s4 | 2 )+V( Am2 " 2 ' Ey =) 2 + 8£;y = Am2 l C/ ^l 2 )' 

U eA = Am U e4 U* 4 y (Am 2 -2£^) 2 +8£;y s Am 2 |c/ S 4| 2 



(19) 



and 



TjM 

U fJ,4: 



2 t / l+4Am 2 |C/ 3 4| 2 (l-|y s 4| 2 )/(2By 3 -Ar7i 2 (l-2|a s4 | 2 )+^/^m 2 -2£;V ii ) 2 +8£;V;Am 2 |;7 s 4| 2 ) 2 

Am U^U* si \J (Am^-2EV a ) 2 +8EV 3 Arn 2 \U s4 .\ 2 ' 

(20) 

The oscillation probability is then derived just as in the standard (3+1) neutrino oscillation scenario (from Eq. [T8|). 



-> v.) = 4\U%\ 2 \U™\ 2 sin 2 (1.27 Am 2 M L/E) , (21) 



where one has replaced Am 2 with Am 2 M from Eq. [T5] and £/ a4 with from Eqs. Q1J] and The effective mixing 
amplitude in matter, sin 2 29ff e = 4|[/ £ Jf| 2 |£/* 4 1 2 , can be expressed in terms of \U e i\ and \U^\ as 



«in 2 9/9 M - 16(Am 2 ) 4 |[/ e4 | 2 |l/„ 4 | 2 |t/ 3 4| 4 

sm — 7 \ 5 

F {(Am 2 -2EVs) 2 +8EVsAm 2 \U B 4\ 2 )(2EV s ~Am 2 (l-2\Us4\ 2 ) + y / (Am 2 -2EV s ) 2 +8EV s Am 2 \Us4\ 2 ) 



(22) 
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Note that sin 2 28ff e reduces to sin 2 28^ = 4|[/ e4 | 2 |i7 /i 4| 2 when V s —> 0, as expected. 

As in the case of standard matter effects, resonances occur when the denominator of Eq. [22] is minimized (42[. It 
is clear that, because of the change in sign in V s for neutrinos versus antincutrinos, we generally expect resonances 
to appear at different energies for neutrinos versus antineutrinos, even when the measurements are performed at the 
same L and E, while underlying vacuum oscillation parameters, Am 2 , \U e 4\, \U^4\ and \U S 4\, are assumed to be the 
same for neutrinos and antineutrinos, as expected by CPT conservation. Upon inspection of the denominator in 
Eq. I22[ given the V s = ±A S definition for neutrinos/antincutrinos, one expects resonances to occur in the case of 
antineutrino but not neutrino oscillations. 

Because sin 2 28^f e and AmJ, arc the parameters measured by experiments by fitting to the usual two-neutrino 
appearance probability formula, given in Eq. II 81 it is instructive to examine sin 2 28^f e and Am^ as a function of E, 
A s , sin 28fj, e , and Am 2 . Figures Q] and [2] illustrate the dependence of Am 2 M (Fig. H]) and sin 2 28ff e (Fig. [2j on the 
underlying vacuum oscillation parameters, sin 2 28 Me and Am 2 , for specific values of E and A s . We have purposefully 
picked neutrino energies E close to the LSND, MiniBooNE low energy and MiniBooNE high energy excess mean 
energies. The A s values have been chosen to span the orders of magnitude considered in our fit. The top set of plots 
in each figure illustrates this dependence for the case of neutrinos, while the bottom set illustrates the antineutrino 
case. From the figures, one can see that as A s — » 0, the effective sin 2 28ff e and Am 2 f approach the vacuum sin 2 28^ e 
and Am 2 parameters. However, as A s turns on, they deviate from the vacuum values, and increasingly so with E. 
Resonances can be identified in lighter regions, where sin 2 28^/ e — > 1. 

III. ANALYSIS METHOD 

We assume that the LSND and MiniBooNE excesses can be described by the oscillation probability in Eq.[5TJ and we 
fit for the vacuum oscillation parameters \U e 4\, |£/ M 4| and Am 2 , and A s , on which Am^ and sin 2 28^f e depend. We set 
| C/ S 4 1 2 = 1 — |C/ e 4| 2 — |t/ M 4| 2 . The vacuum oscillation parameters are allowed to vary freely within 10~ 3 < Am 2 < 100 
eV 2 and sin 2 28^ e < 0.01. The sin 2 28^ ie upper bound is limited by requiring \U e 4\ 2 < 0.05 and |£^i4| 2 < 0.05, assuming 
unitarity of the 3x3 matrix at the 5% level. The matter potential is allowed to vary freely within 10~ 13 < A s < 10~ 9 
cV. The range of A s chosen in our fits has been mainly motivated by the assumption that this new potential should 
lead to observable effects at L/E ~ 1 eV 2 , which is supported by the findings in [29| . For comparison, the standard 
matter effect pottential considered in long-baseline neutrino oscillation experiments is \f2GFn e ~ 10~ 13 eV. 

During the fit, the model parameters |£^4|, \U e 4\, Am 2 , and A s are generated and varied according to a Markov 
Chain (30 . [3l[ \ 2 minimization routine. For each variation, the signal predictions for MiniBooNE neutrino mode, 
MiniBooNE antineutrino mode and LSND are calculated using the oscillation probability in Eq . |2"T1 and then compared 
to the observed excesses in the form of a x 2 . 

The LSND and MiniBooNE data sets used in fits presented in this paper are identical to those in [l4j], with 
the exception of the MiniBooNE antineutrino data set, where we use the higher-statistics, updated results from 
the MiniBooNE — > v e search, corresponding to 5.66xl0 20 protons on target (POT) Note that MiniBooNE 
antineutrino data taking is still in progress. The experiment aims to complete its total antineutrino running at the 
end of this spring, with an estimated final antineutrino sample corresponding to ^lOxlO 20 POT. The fits presented 
in this paper should therefore be updated once the higher statistics antineutrino data set from MiniBooNE becomes 
available. 

The MiniBooNE neutrino data set (MBv) is included in the fits in the form of two side-by-side distributions of 
v e and charged-current quasi-clastic (CCQE) candidate events, each as a function of neutrino energy which has 
been reconstructed assuming CCQE scattering. Using both the and v e CCQE reconstructed events in the fit along 
with their correlations significantly reduces some of the systematic errors related to neutrino flux and cross section. 
The full 200-3000 MeV range of v e CCQE data available is used in the fit. The observed event distributions are 
compared to the corresponding Monte Carlo predicted distributions, and a % 2 is calculated using a covariance matrix 
which includes systematic and statistical uncertainties as well as systematic correlations between the predicted v e 
and distributions. During the fit, we vary the v e distribution according to the sterile neutrino vacuum oscillation 
parameters and A s , but we keep the distribution unchanged, despite the possibility of disappearance in the 
MiniBooNE data. We verify that the best-fit model implies disappearance well below the level constrained by 
the MiniBooNE, SciBooNE, and MINOS disappearance analyses |32|-l34|. and assume that the effect of any 
disappearance on the ability of the CCQE sample to constrain the v e CCQE sample is small. The fit method 
follows the details described in [l4j ]. 

The MiniBooNE antineutrino data set (MBP) is included in the fits in the same way as the MBv data set, in the 
form of two side-by-side distributions of v e and P M CCQE candidate events. The full 200-3000 MeV range of v e CCQE 
data is used in the fit. The MBP data fit method also follows the details described in [l4j |. 
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The LSND data set is included in the fits in the form of a v e event distribution from v e p — s- ne + interactions, as 
a function of five positron energy bins between 20 and 60 MeV. We neglect the higher-energy pion-decay-in-flight 
sample and fit only the decay-at-rest sample which makes this a pure P M — > v e search. Because this is a low-statistics 
sample, the LSND \ 2 function is constructed as a log-likelihood function of the observed data, expected signal, and 
expected backgrounds. The fit method follows the details of (hi. l35l. l36j . 

The x 2 returned by each data set is used to construct the total \ 2 , 

Xtotal = Xmbv + Xmbp + Xlsnd > (23) 

which is minimized in the Markov Chain, and then used to extract the best fit and 90% and 99% confidence level 
(CL) allowed oscillation parameter regions. The CL intervals shown in the figures of this paper correspond to the 
standard A% 2 cuts for two (2) degrees of freedom. 

In order to quantify the statistical compatibility between various data sets under a particular hypothesis, we use the 
Parameter Goodness-of-fit (PG) test introduced in [37|]. This test quantifies the level of agreement between various 
data sets by comparing the minimum x 2 obtained by a simultaneous fit to all data sets, Xmin all' *° ^ ne sum °f * ne 
individual minimum x 2 's obtained by a separate fit to each of the data sets, i.e., 

(24) 

i 

where i runs over the data sets considered in the fit yielding Xmin ail- The is obtained from x 2 p G based on the 
number of common underlying fit parameters, ndfp G , using the standard probability distribution function. 



IV. RESULTS 



A. Fit results with no matter effects: A s — 



In this section, the MiniBooNE and LSND results are examined under a (3+1) oscillation hypothesis with no 
matter effects (A s = 0). The results are to be used as a reference in Sees. II V Bl and fVl Note that results presented in 



this section differ from those reported in |14( because of (a) the updated MBD data set being used here and (b) the 
| C/e4 1 2 < 0.05 and |C/ M 4| 2 < 0.05 fit constraints imposed in this analysis. 

A (3+1) fit to oscillations is unable to reconcile the three signatures. The x 2 -probability for the best-fit parameters 
obtained from a simultaneous fit to all three data sets corresponds to 6.8% (x 2 /ndf = 52.89/39). The compatibility 
between all three data sets is found to be 2.3%, using the PG criterion 

PG = prob{ X 2 PGl ndfp G ) = (11.4,4) = 2.3% . (25) 

The source of incompatibility is demonstrated in Fig. [3] While MiniBooNE antincutrino and LSND (antineutrino) 
results yield contours which overlap in regions of high confidence level, the MiniBooNE neutrino results highly exclude 
those regions, and are preferentially shifted to lower (sin 2 26^, Am 2 ) values. The incompatibility is also illustrated in 
Fig.21 where the phase of the MiniBooNE neutrino observed L/E distribution is significantly shifted relative to that 
of the MiniBooNE antineutrino data set or that expected from the LSND best-fit prediction. 



B. Fit results with matter effects: A s =^ 



In this section, the MiniBooNE and LSND results are examined under a (3+1) oscillation hypothesis with matter 
effects (A s + 1 0). The x 2 -probability for the best- fit parameters obtained from a simultaneous fit to all three data sets 
corresponds to 21.6% (x 2 /ndf = 44.54/38). Relative to the (3+1) fit in Sec. EES the best fit x 2 is reduced by 8.35 
units for one extra fit parameter. The compatibility between all three data sets is found to be 17.4%, using the PG 
criterion 

PG = prob(x 2 PG ,ndfp G ) = (9.0, 6) - 17.4% . (26) 

The x 2 probability for each experiment, at the best-fit parameters found by the joint fit, corresponds to 25.0% (MBi^), 
15.3% (MBP) and 12.3% (LSND), supporting the high compatibility reported above. 

The best-fit signal plus background predictions as a function of energy are shown in Fig. (SJ for each of the three data 
sets. For comparison, the (3+1) best-fit predictions are also overlaid. The (3+1) fit with matter effects accommodates 
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Data Set 


Ami, (eV 2 ) 


sm 20 pe 


MB// 


0.75 


0.0016 


MBv 


0.22 


0.08 


LSND 


0.45 


0.012 



TABLE I: Effective best-fit mixing parameters for the data sets considered in the fit. The parameters have been calculated 
using Embv = 600 MeV, Embb = 700 MeV and E LS nd = 40 MeV. A s corresponds to the best-fit value of 2.0xl0~ 10 eV. 



a smaller fraction of the low energy excess observed in MiniBooNE neutrino mode than the (3+1) fit, but describes 
the high energy data significantly better. In antineutrino mode the fit predicts a considerably larger excess than the 
(3+1) fit, across all energies. In the case of LSND, the matter effect fit performs only marginally better than the 
(3+1) fit. 

The allowed parameters at 90% and 99% CL are shown in Fig. [5] The best-fit parameters, indicated on the figure 
in black stars, correspond to 

sin 2 20 Me = 0.010 , 
Am 2 = 0.47 eV 2 , 
As = 0.2 x lO" 10 eV . (27) 

Note that the Am 2 and sin 2 29 ^ e values quoted above correspond to the vacuum oscillation parameters, and are 
consistent with a fourth, mostly sterile neutrino mass cigenstate with Am^ 3> Am 2 2 , Am^, by construction. The 
effective best-fit parameters for each of the three data sets considered in the fit are summarized in Tab. fl] 



V. QUALITATIVE DISCUSSION OF RESULTS 

The (3+1) fit results arc straightforward to interpret. In the case of no matter effects, the single frequency (Am 2 ) 
involved in a two-neutrino oscillation approximation we have employed in these fits is unable to reconcile the Mini- 
BooNE neutrino mode excess with the MiniBooNE and LSND antineutrino excesses, as they show up at different 
L/E (sec Fig. g]). 

Introducing a matter effect potential which flips sign when going from neutrino to antineutrino oscillations, allows for 
modifying the location (Aro^) and amplitude (sin 2 29^f e ) of the oscillation maximum for neutrinos and antineutrinos 
independently, as well as as a function of E (and L). Thus, oscillation probability measurements performed at the 
same L and E can yield different values for the amplitude and location of the oscillation maximum and corresponding 
observed excess depending on whether they are performed using neutrinos or using antineutrinos. Furthermore, 
measurements at different E can yield excesses which point to different amplitude and oscillation maximum location 
even if they are extracted with the same polarity beam. 

Figure [7] is instructive in understanding how this matter effect leads to the distributions shown in Fig. [5j The top 
panel shows the expected — > v e and — > v e oscillation probabilities at MiniBooNE, as a function of neutrino 
energy. The gray line corresponds to the best-fit oscillation parameters but with A s = 0, while the blue solid and 
dashed lines correspond to the best-fit parameters and A s = 2.0 x 10~ 10 eV. Both oscillation probabilities suggest an 
excess at low energy; however, a non-zero A s value allows for oscillations at higher energy at the ^1% level in the 
case of antineutrinos, and lack of oscillations at higher energy in the case of neutrinos, consistend with MiniBooNE 
observations. The bottom panel shows the expected — > v e oscillation probability at LSND. Here, a ~1% oscillation 
probability is preserved with or without a non-zero A s , due to the much lower energy range. 

Also interesting to explore are the ^ and v e disappearance probabilities expected in this oscillation framework, 

P^^v^) = sin 2 26%, sin 2 (1.27 Am 2 M L/E) , (28) 

which can be calculated using 

sin 2 2C=4|^ 4 T(l-|^i| 2 ), (29) 

and the definitions in Eqs. [Inland I2TJ1 In the case of MiniBooNE, the expected and disappearance probabilities 
corresponding to the best-fit parameters of Eq.[57]are shown in Fig. [5] (top). In both cases (neutrino and antineutrino), 
one expects disappearance probabilities on the order of a few percent, on average. The oscillation probability peaks 
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at lower ener gies , where MiniBooNE and P M disappearance searches become dominated by flux and cross section 
uncertainties |32l. l33j . The upcoming joint MiniBooNE/SciBooNE i/„ disappearance search, however, may have some 
sensitivity to this effect. Similarly, the bottom panel of Fig. [8] shows the v e disappearance probability expected at 
LSND, corresponding to the best- fit parameters of Eq. [57] The expected oscillation probability for E ~ 30 — 40 MeV 
is approximately 10%. The KARMEN and LSND ^ e -Carbon charged-current cross section measurements which are 
available in this energy range do not have sufficient sensitivity to address this j39l | . 

Figure |H] is instructive in understanding how the effective oscillation parameters vary with neutrino energy. The 
top panel of Fig. [S] shows the effective Am^ oscillation parameter in matter as a function of neutrino energy, for the 
best-fit oscillation and A s values. The solid black line indicates the Am^ f dependence on E in the case of a neutrino 
beam; the dashed line shows the antineutrino dependence. The mean neutrino energy for each data set is also shown 
by a vertical line. Similarly, the bottom panel of Fig. [§] shows the effective sin 2 28jf e oscillation parameter in matter, 
and how it varies with neutrino energy, assuming the best-fit parameters obtained in Sec. IIVBI It is evident from this 
figure that, in the case of neutrino oscillation experiments performed in the few MeV to few GeV range, one should 
not expect large appearance amplitudes, since the effective mixing amplitude is always less than that expected from 
the underlying vacuum parameter, sin 2 29^ e = 0.010. On the other hand, one expects to observe resonance-like effects 
with antineutrino experiments performed in the > 1 GeV range. 

As a final discussion point, we remark on the size of A s . While the form of the matter potential we have assumed 
is the same as that of the standard matter effect in three-neutrino oscillations (Asjvf=constant), the resulting best-fit 
value suggested by the MiniBooNE and LSND data sets is roughly 500 times larger than that of standard matter 
effect. The result is consistent with the analytical picture presented in (2j|, where the authors state that in order to 
obtain observable, non-negligible matter effects in a model with a single sterile neutrino, one needs a matter effect 
potential of A s ~ 10~ 10 cV. The large value of A s may be quite difficult to explain in simple extensions to the theory. 
However, the good agreement of all three signatures within this three-parameter model motivates the development of 
theoretical interpretations which would lead to such a phenomenological effect. 



VI. EXPERIMENTAL IMPLICATIONS FOR SOLAR AND ATMOSPHERIC OSCILLATIONS 

Even though the vacuum oscillation parameters allowed in this model are constructed so as not to interfere with 
the atmospheric and solar oscillation scales, the E and V s dependence can drive the effective Am 2 M and sin 2 26^fg 
parameters to the solar and atmospheric parameter regions. Therefore, one must ask whether constraints arise from 
solar/reactor and atmospheric/long-baseline accelerator experiments. 

Reactor experiments pose little to no constraints to this scenario, as the low energy assures that the effective 
oscillation parameters deviate very little from the vacuum parameters (which, according to the best fit, correspond 
to Am 2 = 0.47 eV 2 and sin 2 28 ee = 0.19). In fact, the reactor anomaly recently identified in short-baseline reactor 
experiments is consistent with small-amplitude oscillations due to a heavy (>1 eV 2 ) mostly-sterile mass eigenstate 
[ll| . and could reasonably be accommodated within the matter effect scenario. Joint fits including the reactor data 
are necessary for more quantitative tests. 

It is likely that such matter effects affect neutrino propagation in the sun. However, this is probably dependent on 
the nature of V s . Such effects should be investigated for particular V s underlying scenarios. 

In the case of atmospheric/long-baseline accelerator oscillation experiments, we use the MINOS experiment as an 
example to investigate the effects of the best-fit model in Sec. IIVBI MINOS employs a near/far detector setup at 
~ 0.75/735 km from a neutrino source, in order to look for and disappearance driven by the atmospheric Am 2 
scale, Am 2 2 - The MINOS neutrino energy spectrum spans the 1-10 GeV energy range, with a peak neutrino energy at 
3 GeV. Because the Am^ deviation from the vacuum oscillation parameter becomes larger, overall, with increasing 
energy, and MINOS sits at a much higher energy than MiniBooNE and LSND, one might expect noticeable effects. 
From Fig. [§] one sees that for E ~ 3 GeV, the effective Am\ [ for neutrinos is similar to or larger than the vacuum 
Am 2 , and therefore one might expect any possible oscillations to average out to half the disappearance amplitude, 
sin 2 26^ /2, for both the near and far MINOS detectors. Here, sin 2 29^ corresponds to 

S in 2 2^ = 4|f/^| 2 (l-|^| 2 ), (30) 

which is the standard disappearance probability amplitude for (3+1) sterile neutrino oscillations, except with 
— > UM, from Eq. [20] Figure [TOl shows the sin 29^ dependence on E for the best-fit parameters. One sees 
that sin 2 29ff^ is always equal to or smaller than the vacuum oscillation parameter sin 2 29^ = 0.19 for any neutrino 
oscillation experiment. Therefore, the MINOS neutrino mode data set is insensitive to this model, since systematic 
uncertainties on the absolute event rate prediction at the near detector arc much larger than this level. 

In the case of antineutrino oscillations, however, the Am^ value decreases with E and becomes negative beyond 
Eq = Am 2 / (2A S ) ~ 1180 MeV. Note, however, that sin 2 (1.27Am 2 / L/i?) is insensitive to the sign of Am 2 M . Beyond 
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this point, therefore, the effective |Am 2 /f | starts increasing linearly with E. This is illustrated in Fig. [9l Given 
the best-fit oscillation parameters, we see that (Am 2 ^) becomes close to the atmospheric Am^ for a narrow range 
of energies around EL 1200 McV. This range is on the very low edge of the MINOS energy range, and while the 
sin 2 28^ value for antineutrinos at this energy is predicted to be nearly maximal, without an actual fit it is difficult 
to make quantitative statements as to the level of constraint provided by MINOS antineutrino data, since the Am^ 
and sin 2 28jf^ vary rapidly in the MINOS antineutrino energy range. On the other hand, it is reasonable to expect 
differences in the neutrino and antineutrino oscillation parameters obtained from MINOS and Un disappearance 

searches, respectively, such as those presented in [I(|[4l|. The bottom panels of Fig. ITTI show the observable v ^ 
disappearance probabilities expected at MINOS in neutrino and antineutrino mode at both the far and near detectors. 

The top panels of Fig. lTTI show the corresponding observable v ^— > v e appearance probabilities expected at MINOS 
near and far detectors [43|. In neutrino mode, one expects a < 1% appearance oscillation probability across the 1- 
10 GeV energy range at the far detector, which is far beyond MINOS' sensitivity, but in antineutrino mode, on expects 
large appearance effects which may be measurable in a far to near comparison, given enough statistics. 



VII. CONCLUSIONS 



We have examined MiniBooNE and LSND results in oscillation fits to a model with a single, mostly-sterile neutrino 
mass eigenstate at a Am 2 ~ 1 eV 2 , with and without the presence of an effective matter- like potential of the form 

V s = ±A S , (31) 

experienced only by sterile neutrino/antineutrino states. We find that the compatibility among LSND and MiniBooNE 
neutrino and antineutrino data sets increases significantly in the presence of a non-zero A s , with the best- fit parameters 
corresponding to the vacuum mixing parameters sin 2 28^ = 0.010 and Am 2 = 0.47 cV 2 , and A s = 2.0 x 10~ 10 cV. The 
best-fit parameters are consistent with reactor long-baseline, atmospheric, and accelerator long-baseline (neutrino) 
results, and can reasonably accommodate the recent reactor short-baseline anomalous result. Implications for the 
MINOS antineutrino data set have also been considered; the best-fit model predicts some observable effects in the 
MINOS antineutrino samples, and therefore MINOS' sensitivity to those effects should be explored further. We invite 
phenomcnological interpretations of this model. 
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FIG. 1: Effective Amjj (colored contours, with color scale on the right) as a function of the vacuum oscillation parameters Am 2 
and sin 2 28 Me for different sets of neutrino energy values E and matter effect potential A s . The top set of plots corresponds to 
neutrinos. The bottom set corresponds to antineutrinos. 
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Am 2 and sin 2 26^ for different sets of neutrino energy values E and matter effect potential A s . The top set of plots corresponds 
to neutrinos. The bottom set corresponds to antineutrinos. 
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FIG. 4: The MiniBooNE neutrino (top) and antineutrino (bottom) observed excess distributions each as a function of L/E, 
E being the reconstructed neutrino energy, E® E . The error bars include systematic and statistical uncertainties. Overlaid are 
signal predictions corresponding to the LSND (3+1) best-fit parameters. The figure is taken from [2r| . 
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overlaid on the MiniBooNE neutrino (top), antineutrino (middle), and LSND (bottom) observed spectra. The no oscillations 
prediction is also shown, in gray. In the case of the MiniBooNE neutrino distributions, the matter effect fit predicts a lower 
excess at low energy, but also no excess at higher energy, which reduces the tension otherwise present in the (3+1) fit. In the 
case of the MiniBooNE antineutrino distributions, the matter effect fit predicts both low and high energy excess. In the case 
of the LSND distributions, the two fits are essentially indistinguishable. 



15 



□ 90% CL 

□ 99%CL 



sirrse 



lie 



pi— 
> 

<u 10 

CM 

E 

< 



10 1 
10 



□ 90% CL 
n99%CL 



TO' 12 



10"" 



10" 



A s {eV) 



FIG. 6: The 90% and 99% confidence level (CL) allowed vaccum oscillation parameters and A s , obtained from a joint (3+1) 
with matter effect fit to MBv, MBi? and LSND data sets. 



16 




FIG. 7: Effective appearance oscillation probabilities for MiniBooNE (top) and LSND (bottom) as a function of neutrino 
energy, calculated using the best-fit values in Eq. 1271 The gray lines correspond to the same oscillation parameters but with 
A a set to zero. 
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FIG. 8: Effective disappearance oscillation probabilities for MiniBooNE (top) and LSND (bottom) as a function of neutrino 
energy, calculated using the best-fit values in Eq. 1271 The gray lines correspond to the same oscillation parameters but with 
A a set to zero. 
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FIG. 9: The effective Am\j and sin 2 26^, mixing parameters corresponding to the best-fit values in Eg. 1271 each as a function of 
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appearance (top) and v M disappearance (bottom) oscillation probabilities for MINOS as a 
function of neutrino energy, calculated using the best-fit values in Eq. 1271 The left figures show neutrino mode oscillation 
probabilities. The right figures show antineutrino mode oscillation probabilities. The gray lines correspond to the same 
oscillation parameters but with A 3 set to zero. 



